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Abstract: Very high energy (VHE; > 100 GeV) 7-ray observations of PG 1553+1 13 were made with the 
High Energy Stereoscopic System (HESS) in 2005 and 2006. A strong signal, ~10 standard deviations, 
is detected by HESS during the 2 years of observations (24.8 hours live time). The time-averaged energy 
spectrum, measured between 225 GeV to ~ 1.3 TeV, is characterized by a very steep power law (photon 
index of V = 4.5 ± 0.3 stat ± 0.1 syst ). The integral flux above 300 GeV is -3.4% of the Crab Nebula 
flux and shows no evidence for any variations, on any time scale. H+K (1.45— 2.45^im) spectroscopy of 
PG 1553+113 was performed in March 2006 with SINFONI, an integral field spectrometer of the ESO 
Very Large Telescope (VLT) in Chile. The redshift of PG 1553+1 13 is still unknown, as no absorption or 
emission lines were found. 



Introduction 

Evidence for VHE (>100 GeV) 7-ray emission 
from the high-frequency -peaked BLLac object 
PG 1553+113 was first reported by the HESS 
collaboration [3] based on observations made in 
2005. This detection was later confirmed [6] with 
MAGIC observations in 2005 and 2006. The 
measured VHE spectra are unusually soft (pho- 
ton index T=4.0±0.6 and T=4.2±0.3 for the HESS 
and MAGIC experiments, respectively) but the er- 
rors are large, clearly requiring improved measure- 
ments before detailed interpretation of the com- 
plete SED is possible. Further complicating any 
SED interpretation is the absorption of VHE pho- 
tons [10, 11] by pair-production on the Extra- 
galactic Background Light (EBL). This absorption, 
which is energy dependent and increases strongly 
with redshift, distorts the VHE energy spectra ob- 
served from distant objects. For a given redshift, 
the effects of the EBL on the observed spectrum 
can be reasonably accounted for during SED mod- 
eling. Unfortunately, the redshift of PG 1553+1 13 
is unknown, despite many attempts to measure it 
(see, e.g., [8, 16]). 



In 2005 and 2006, a total of 30.3 hours of HESS 
observations were taken on PG 1553+113. The 
2005 HESS observations are exactly the same as 
reported in [3]. The good-quality exposure is 24.8 
hours live time. The data are processed using the 
standard HESS calibration [2] and analysis tools 
[7]. Soft cuts [3] are applied to select candidate 
7-ray events, resulting in an average post-analysis 
energy threshold of 300 GeV at the mean zenith 
angle of the observations, 37°. 

HESS Results 

A significant VHE 7-ray signal is detected in each 
year of HESS data taking. The total observed ex- 
cess is 785 events, corresponding to a statistical 
significance of 10.2 standard deviations (a). Ta- 
ble 1 shows the results of the HESS observations of 
PG 1553+113. Figure 1 shows the on-source and 
normalized off-source distributions of the square 
of the angular difference between the reconstructed 
shower position and the source position (6 2 ) for 
all observations. The background is approximately 
flat in 2 as expected, and there is a clear point- 
like excess of on-source events at small values of 
9 2 , corresponding to the observed signal. The peak 



Multi-wavelength Observations of PG 1553+1 13 




Figure 1: The distribution of 6 for on-source 
events (points) and normalized off-source events 
(shaded) from observations of PG 1553+1 13. The 
dashed curve represents the 6 2 distribution ex- 
pected for a point source of VHE 7-rays at 40° 
zenith angle with a photon index T = 4.5. The 
vertical line represents the cut on 8 2 applied to the 
data. 

of a two-dimensional Gaussian fit to a sky map of 
the observed excess is coincident with the position 
of PG 1553+113. 

The photon spectrum for the entire data set is 
shown in Figure 2. These data are well fit (x 2 
of 8.4 for 5 degrees of freedom) by a power-law 
function, &NIA.E ~ E~ r , with a photon index 
r = 4.5 ± 0.3 s tat ± 0.1 S y S t. Fits of either a power 
law with an exponential cut-off or a broken power 
law do not yield significantly better x 2 values. 

The observed integral flux above 300 GeV for the 
entire data set is I(>300 GeV) = (4.6 ± 0.6 sta t ± 



Table 1 : Shown are the excess, the significance of 
the excess, and the integral flux above 300 GeV, 
from HESS observations of PG 1553+113. The 
flux units are 10~ 12 cm~ 2 s . The systematic er- 
ror on the flux is 20% and is not shown. 



Epoch Time Excess Sig I(>300GeV) 

[h] [a] [f.u.] 

2005 7.6 249 6.0 5.44 ±1.23 

2006 17.2 536 8.3 4.22 ± 0.72 



Total 24.8 785 10.2 4.56 ± 0.62 



Figure 2: The overall VHE energy spectrum ob- 
served from PG 1553+1 13. The dashed line repre- 
sents the best x 2 nt of a power law to the observed 
data. The upper limits are at the 99% confidence 
level [9]. Only the statistical errors are shown. 

0.9 sys t) x 10~ 12 cm _2 s _1 . This corresponds 
to -3.4% of I(>300 GeV) determined from the 
HESS Crab Nebula spectrum [5]. The integral flux, 
I(>300 GeV), is shown in Table 1 for each year of 
observations. Figure 3 shows the flux measured for 
each dark period. There are no indications for flux 
variability on any time scale within the HESS data. 
The data previously published [3] for HESS ob- 
servations of PG 1553+1 13 in 2005 were not cor- 
rected for long-term changes in the optical sensitiv- 
ity of the instrument. Relative to a virgin telescope, 
the total optical throughput was decreased by 29% 
in 2005 and 33% in 2006. Correcting [5] for this 
decrease, using efficiencies determined from simu- 
lated and observed muon events, increases the flux 
measured from the object. The effect of this cor- 
rection is larger for soft spectrum sources than it is 
for hard spectrum sources. Due to the correction, 
the flux measured in 2005 is three times higher 
than previously published. 

On July 24 and 25, 2006, PG 1553+113 
was observed by the Suzaku X-ray satel- 
lite ([14], Suzaku Observation Log: 
http://www.astro.isas.ac.jp/suzaku/index.html.en). 
On these two dates HESS observed PG 1553+113 
for 3.1 hours live time, resulting in a marginally 
significant excess of 101 events (3.9cr). The 
average flux measured on these two nights is 
I(>300 GeV) = (5.8 ± 1.7 stat ± 1.2 syst ) x 10~ 12 

— 9 — 1 

cm s . 
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The H+K-band spectrum of PG 1553+113 mea- 
sured here is shown in Figure 4. The signal-to- 
noise reach is ^250 in the H-band and ^70 in the 
K-band. The observed near-IR spectrum is fea- 
tureless apart from some residuals from the atmo- 
spheric corrections. Thus in neither the broad-band 
images nor in the spectrum are the influences of 
the gas of a host galaxy or the AGN detected, even 
though PG 1553+113 is bright in the IR. As a re- 
sult, a redshift determination from these observa- 
tions is not possible. 



Figure 3: The integral flux, I(>300 GeV), mea- 
sured by HESS from PG 1553+113 during each 
dark period of observations. The horizontal line 
represents the average flux for all the HESS obser- 
vations. Only the statistical errors are shown. 
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Figure 4: The measured H+K-band spectrum of 
PG 1553+1 13 in relative flux units. The gap is due 
to the highly reduced atmospheric transmission be- 
tween H and K bands. 



Discussion 

As the redshift of PG 1553+113 is likely z > 0.2, 
the observed VHE spectrum should be strongly af- 
fected by VHE 7-ray absorption on the EBL. If the 
redshift were known the spectrum intrinsic to the 
source could be reconstructed assuming a model 
of the EBL density. However, the EBL SED is not 
well-determined. Using a Maximal EBL model at 
the level of the upper limits from [4] or a Mini- 
mal model near the EBL lower limits from galaxy 
counts [12] can yield a significantly different in- 
trinsic spectrum. Figure 5 shows the intrinsic spec- 
trum versus redshift for both the Maximal and Min- 
imal EBL parameterizations. Here, scaled models 
of [13] are used, exactly as described in [3]. The 
redshift of the AGN can be limited using assump- 
tions for the intrinsic spectrum. Assuming the in- 
trinsic photon index is not harder than Tint = 1.5, 
a limit of z < 0.69 is thus determined from the 
Minimal EBL model. 



SINFONI Near-IR Spectroscopy 

The determination of the redshift of an AGN is 
generally based upon the detection of emission or 
absorption lines in its spectrum. In an attempt to 
detect absorption features from the host galaxy or 
emission lines from the AGN, H+K (1 .50-2.40yum) 
spectroscopy of PG 1553+1 13 was performed with 
SINFONI, an integral field spectrometer mounted 
at Yepun, Unit Telescope 4 of the ESO Very Large 
Telescope in Chile. The source was observed on 
March 9, 2006 and March 15, 2006. The resulting 
images are spatially unresolved and no underlying 
host galaxy is detected. 



Conclusion 

With a ^25 h data set, ^3 times larger than 
previously published [3], the HESS signal from 
PG 1553+113 is now highly significant (<~10<r). 
Thus the evidence for VHE emission from 
PG 1553+113 previously reported is now clearly 
verified. However, the flux observed in 2005 is 
now ^3 times higher than initially reported due 
to an improved calibration of the absolute energy 
scale of HESS, and agrees well with the flux mea- 
sured in 2005 by MAGIC [6]. The statistical er- 
ror on the VHE photon index is still rather large 
(~0.3), primarily due to the extreme softness of 
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Figure 5: The photon index r; nt determined from 
a power-law fit to the intrinsic spectrum of PG 
1553+1 13 (i.e. the H.E.S.S. data de-absorbed with 
an EBL model) for a range of redshifts. The con- 
tours reflect the la statistical uncertainty of the fits. 
The upper curve is the sum of r int for the Mini- 
mal" model and twice its statistical error. 

the observed spectrum (T = 4.5). The total HESS 
exposure on PG 1553+1 13 is ^25 hours. Barring a 
flaring episode, not yet seen in two years of obser- 
vations, a considerably larger total exposure (~100 
hours) would be required to significantly improve 
the spectral measurement. However, the VHE flux 
from other AGN is known to vary dramatically 
and even a factor of a few would reduce the ob- 
servation requirement considerably. Should such 
a VHE flare occur, not only will the error on the 
measured VHE spectrum be smaller, but the mea- 
sured photon index may also be harder (see, e.g., 
[1]). Both effects would dramatically improve the 
redshift constraints and correspondingly the accu- 
racy of the source modeling. Therefore, the VHE 
flux from PG 1553+1 13 will continue to be moni- 
tored by HESS. In addition, the soft VHE spectrum 
makes it an ideal target for the lower-threshold 
HESS Phase-II [15] which should make its first ob- 
servations in 2009. 
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